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^SJ ' Constrained Relativistic Mean Field (RMF) calculations have been 

carried out to estimate excitation energies relative to the ground state 

Qv i for superdeformed bands in the mass regions A ~ 190 and A ~ 150. 

It is shown that RMF theory is able to successfully reproduce the 
recently measured superdeformed minima in Hg and Pb nuclei. 

CJ I Superdeforination has become the recent years one of the most interest- 

^ I ing topics of nuclear structure studies. After the well known area in the 

vicinity of the mass number A ~ 150 a second region with A ~ 190 has 
been discovered, where an impressive number of results has been obtained. 
However, despite the rather large amount of experimental information on 
5^ I superdeformed bands, there are still a number of very interesting proper- 

ties, which have not yet been measured. A characteristic example is the 
excitation energy between ground and superdeformed bands. The excitation 
energy and the well depth of the superdeformed minimum are amongst the 
most important factors which affect the decay of the superdeformed bands 
to the ground state. 

Very recently discrete 7 rays directly connecting states of a superdeformed 
band in ^^''Hg to the yrast states have been discovered |jl|. This has made 
it possible to determine accurately the excitation energy of the levels of the 
superdeformed band above the yrast line. Extrapolating to zero angular mo- 
mentum the superdeformed minimum was found to be 6.017 MeV above the 
ground state. Similar measurements at about the same time have also been 
reported for the superdeformed band of ^^"^Pb nucleus [|] . The excitation en- 
ergy between the SD band and the ground state was estimated to be 4.471(6) 
MeV. 



There are several non-relativistic theoretical predictions for the excitation 
energy relative to the ground state. Such as Hartree-Fock calculations with 
density dependent Skyrme [Q or Gogny interactions ^ as well as calculations 
using the Strutinsky method built on a Woods-Saxon average potential [|]. 
In this letter we report and discuss the predictions of the RMF theory for 
the excitation energies relative to the ground state of the SD bands of ^^^Hg 
and ^^^Pb. In addition predictions for nuclei in the rare earth mass region 
are also provided. 

Relativistic Mean Field (RMF) theory has recently gained considerable 
success in describing various facets of nuclear structure properties. With a 
very limited number of parameters, RMF theory is able to give a quantitative 
description of ground state properties of spherical and deformed nuclei p, Q 
at and away from the stability line [l^, |ll|, |l^. Moreover good agreement 
with experimental data has been found recently for collective excitations such 



as giant resonances [^ and for twin bands in rotating superdeformed nuclei 



T^, ITB]. Afansajev et al, have recently, and with great success carried out a 



systematic investigation of the entire A~140 to A~150 mass region [^ . 
The starting point of RMF theory is a standard Lagrangian density |^ 

£ = ^ ('j(id - g^u - gppr - eA) - m - g„a) tp 

-^R.^R'^'^ + ^m^-^-lF^.F'^^ (1) 

which contains nucleons i/j with mass m, a-, u-, p-mesons, the electromag- 
netic field and non-linear self-interactions of the cr-field, 

Uia) = \mla' + \g,a' + \g,a\ (2) 

In the present work we have carried out constrained RMF calculations 



by imposing a quadratic constraint H' [|T8|, for the quadrupole moment us- 
ing a mean field Hamiltonian Yi^MF and minimizing instead of Yi^MF the 
expectation value of YYrmf + H' where 

H' = |(< Q > -q? (3) 

and q is the actual value of the quadrupole moment. This constraint enables 
us to trace the energy surface as function of the quadrupole moment. 

For our calculations we have adopted the frequently used parameter set 

us- 



NLl |T^, ^, along with the recently proposed parametrization NL3 pi 
ing a new version of the "axially deformed" code [^]. The values of the two 
parameter sets are listed in Table 1. In the calculations pairing has been 
included within the BCS formalism. The experimental odd-even mass differ- 
ences were used to determine for the ground state gap parameters Ap (for 



protons) and A„ (for neutrons). These values were then used to determine 
the corresponding monopole pairing constants Gp and Gn using the energy 
cut-off of 82^4"^/^ MeV. In order to treat properly the pairing correlations 
on going to the super deformed minimum the corresponding gaps were not 
kept constant but calculated using these strength parameters for each value 
of the constraint. 

In Fig. 1 we show the potential energy surface for the nucleus ^^^Hg as a 
function of the quadrupole parameter q in Eq. (^ , calculated with the param- 
eter sets NL3 and NLl. It is seen that for both Lagrangian parametrizations 
the RMF theory predicts an oblate shape for the the ground state of the 
nucleus ^^^Hg . This is in accordance with experiment p3[. The calculations 



predict in addition a stationary point at prolate deformation. In an axially 
symmetric calculation it looks like a local minimum, but including triax- 
ial deformations it corresponds probably to a saddle point connected to the 
oblate minimum by a valley in the {3-^ plane, which avoids the maximum at 
/3 = 0. For this nucleus the excitation energy E^ of the superdeformed mini- 
mum with respect to the oblate ground state is 5.62 MeV for the parameter 
set NLl. It is close to the measured value (6.02 MeV) |l|. The parameter set 
NL3 does even better giving a value of 6 MeV in excellent agreement with the 
experiment. The depth of the potential well of the isomeric superdeformed 
state is about 2.66 MeV for NLl and 1.25 MeV for NL3. 

We have also performed similar calculations for the nucleus ^^^Hg. These 
are shown in Fig. 2 . For this nucleus the excitation energy E^ is 3.93 MeV 
for NLl while 4.46 MeV for NL3. The depths of the isomeric wells were 
found to be 2.70 MeV and 1.17 MeV respectively. 

It is seen from Figs. 1 and 2 that going from ^^^Hg to ^^^Hg the excitation 
energy E^ as well as the depth of the superdeformed minimum are decreased 
for both parametrizations. For NLl the relative decrease is 30 % and 1.5 % 
respectively while for NL3 it is 26 % and 6.5 %. Moreover, it is also seen that 
for the two Hg nuclei NL3 predicts higher excitation energies and shallower 
potential wells as compared with NLl. An interesting feature which is com- 
mon for both parametrizations is the fact that the product of the well depth 
and width of the superdeformed barrier decreases by about 20% on going 
from ^^^Hg to ^^^Hg. This quantity is relevant for the tunneling probability 
between the superdeformed and the other minima and therefore for the life 
time in the superdeformed well. Similar observations have also been made 
in experimental investigations [^]. It should be noted at this point that the 
parameter set NL3 gives values for the ground state binding energies which 
are in excellent agreement with the empirical values. Specifically NL3 pre- 
dicts for i^^Hg 1519.87 MeV and for ^^^Hg 1535.90 MeV the experimental 
values being 1519.20 MeV and 1535.50 MeV respectively. The corresponding 
predictions of NLl are 1524.99 and 1540.31 MeV respectively. 

Non-relativistic Density Dependent Hartree-Fock calculations predict val- 



ues for the excitation energies E^ of ^^^Hg, which deviate somewhat from the 
experimental values. For Skyrme forces one finds 0] 5.0 MeV, while the 
Gogny's force predicts 6.9 MeV. Finally, using a Woods-Saxon potential 
[^] a value of 4.6 MeV was obtained. We can conclude that the predictions of 
RMF theory are in better agreement with experiment than other theoretical 
studies. 

Our calculations stayed strictly in the mean field approximation. There 
were no corrections of spurious rotational contributions subtracted. Those 
corrections would influence in some way or another the energy of all the 
stationary points in the energy surface. However only the differences between 
the minima would be of importance. To calculate these contributions in an 
appropriate way is a very difficult task, which definitely goes beyond the 
present state of the art. 

We next, carried out calculations for the nucleus ^^^Pb using the NL3 
force. The potential energy landscape is shown in Fig. 3 (left side). RMF 
theory predicts for the excitation energy E^ a value 4.53 MeV, which is very 
close to the experimentally measured value 4.471(6) MeV. The depth of the 
isomeric well was found to be 2.2 MeV. The Hartree-Fock calculations with 
a Skyrme interaction give an estimate of 4.86 MeV, which is also in good 
agreement with the measured value. The Woods-Saxon potential gives a 
lower value of 3.8 MeV. In contrast to the Hg-imdei we observe here a rather 
sharp maximum between the ground state minimum at zero deformation and 
the superdeformed minimum. However, this does not mean pairing correla- 
tions vanish at this saddle point. As expected we observe a local increase 
of pairing in going over the saddle because of the increasing level density at 
this point. 

Recently its was observed |^ that the superdeformed band of ^^^Pb is 
populated at lower spin values than that of ^^^Hg and de-excites towards 
normally deformed states at a spin value lower than that of ^^^Hg. This fact 
suggests that the barrier separating the superdeformed and normal deformed 
minima in the potential energy surface is higher in ^^^Pb than in ^^^Hg. This 
agrees with the predictions of the RMF theory. Our calculations using the 
force NL3 for these two isotonic nuclei (N=90) show that the height of the 
barrier of the nucleus ^^^Pb is about 1 MeV higher than ^^^Hg and also the 
excitation energy E^ of ^^"^Pb is slightly higher than that of ^^^Hg. 

Superdeformed nuclei in the rare earth region are of great interest and 



there are projects for the measurement of the excitation energies E^ [^ 
We have therefore, also performed calculations for some even-even rare earth 
nuclei which exhibit superdeformed bands. In Table 2 we list the predictions 
of RMF theory for the excitation energy E^. It is seen for Gd and Dy isotopes 
the Ex values increase with increase of the mass number. However, the E^ 
energies of ^'^^Gd and ^^^Dy nuclei appear to be larger than their neighboring 
counterparts. This could attributed to the magic character of these nuclei 



which are expected to have a deeper minimuni. Moreover, it is observed 
that the isotonic nuclei, ^^^Sm, ^^^Gd, (N=80) , ^^^Gd, ^^"Dy (N=84), ^^^Gd, 
^^^Dy (N=86) have rather similar excitation energies E^. 

In the same table we show predictions of RMF theory for the ground state 
(g.s) binding energy (BE) together with the corresponding empirical values 
(BEea-p ), taken from the most recent compilation of Audi and Wapstra [^ 



It is seen that RMF theory with the effective force NL3 is able to reproduce 
the experimental values with very high accuracy, the disagreement being less 
than 0.1 %. The excellent predictions for the g.s binding energies give us 
confidence for the correctness of our estimate for the excitation energy iJ^,. 
Finally, we show in Fig. 3 (right side) as an illustration the energy curve 
of the nucleus ^"^^Gd calculated with the parameter set NL3. The excitation 
energy E^ is 10.63 MeV and the depth of the superdeformed minimum is 
0.66 MeV. 

In summary, we have performed constrained RMF calculations using the 
effective forces NLl and NL3 for several nuclei in the A ~ 190 and A ~ 
150 regions to estimate the excitation energy of the superdeformed minimum 
relative to the ground state. It has been shown that RMF theory is able to 
reproduce the recently reported values for the excitation energies of ^^^Hg 
and ^^"^Pb nuclei with a higher precision than other non-relativistic calcula- 
tions. Moreover predictions for the heights of the superdeformed barrier in 
the A ~ 190 mass region have also been given. The calculated values agree 
with scenarios suggested by experimental finding. Finally, RMF theory pre- 
dictions have been also made for the excitation energy E^ of the well known 
superdeformed nuclei in the rare earth mass region, which should also be 
useful for comparisons with future experiments. 

This work has been supported by the Bundesministerium fiir Bildung und 
Forschung under the project 06 TM 785. We thank R. M. Lieder for useful 
discussions. 
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Table 1: The parameter sets of the effective interactions NLf and NL3 of the 
Lagrangian of the RMF theory used in the present work. The masses are 
given in (MeV) and the couphng constant g2 in (fm~^) 



M m„ m^ nip g„ g^ g^ g2 ^3 

NLl 938 492.250 795.359 763.0 10.138 13.285 4.976 -12.172 -36.265 
NL3 939 508.194 782.501 763.0 10.217 12.868 4.474 -10.431 -28.885 



Table 2: The predictions of the RMF theory for the excitation energy E^ 
relative to the ground state of some Sm, Gd and Dy isotopes using the 
effective force NL3. Also shown are the calculated and empirical values of 
the (g.s) binding energies. 



Nucleus 


BE 


BEe^rp 


E, 


"^Sm 


1177.02 


1176.62 


6.80 


i^^Gd 


1184.62 


1184.12 


6.45 


i^^Gd 


1205.43 


1204.44 


10.63 


i^^Gd 


1220.91 


1220.77 


7.82 


isoGd 


1236.85 


1236.40 


8.00 


i48Dy 


1211.47 


1210.75 


10.32 


150 J3y 


1228.35 


1228.39 


7.98 


152Dy 


1245.69 


1245.33 


8.32 


l54Dy 


1262.94 


1261.75 


9.90 



Figure Captions 

Fig. 1 Potential energy surface as a function of the quadrupole parameter 
q for the nucleus ^^^Hg calculated using two different parameter sets 
NL3 (upper part) and NLl (lower part). 

Fig. 2 Potential energy surface as a function of the quadrupole parameter q 
for the nucleus ^^^Hg calculated with two different parameter sets NL3 
(upper part) and NLl (lower part). 

Fig. 3 Potential energy surface as a function of the quadrupole parameter 
q for the nuclei ^^^Pb (l.h.s) and ^^^Gd (r.h.s.) 
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Gd E = 10.63 MeV V= 0.66 MeV 
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